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ABSTRACT 

We report spectroscopic and interferometric observations of the moderately 
metal-poor double-lined binary system HD 195987, with an orbital period of 
57.3 days. By combining our radial-velocity and visibility measurements we de- 
termine the orbital elements and derive absolute masses for the components 
of Ma = 0.844 ± 0.018 Mq and Mb = 0.6650 ± 0.0079 Mq, with relative 
errors of 2% and 1%, respectively. We also determine the orbital parallax, 
TTorb = 46.08 ±0.27 mas, corresponding to a distance of 21.70 ±0.13 pc. The par- 
allax and the measured brightness difference between the stars in V, H, and K 
yield the component absolute magnitudes in those bands. We also estimate the 
effective temperatures of the stars as T^^ = 5200 ±100 K and T^f = 4200 ±200 K. 
Together with detailed chemical abundance analyses from the literature giving 
[Fe/H]?5i —0.5 (corrected for binarity) and [a;/Fe]= ±0.36, we use these phys- 
ical properties to test current models of stellar evolution for metal-poor stars. 
Among the four that we considered, we find that no single model fits all observed 
properties at the measured composition, although we identify the assumptions 
in each one that account for the discrepancy and we conclude that a model with 
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the proper combination of assumptions should be able to reproduce all the ra- 
diative properties. The indications from the isochrone fits and the pattern of 
enhancement of the metals in HD 195987 are consistent with this being a thick 
disk object, with an age of 10-12 Gyr. 

Subject headings: binaries: spectroscopic — stars: fundamental parameters — 
stars: abundances — stars: individual (HD 195987) 

1. Introduction 

Accurate determinations of the physical properties of stars in binary systems (mass, 
radius, temperature, etc.) provide for fundamental tests of models of stellar structure and 
stellar evolution. The most basic of those stellar properties is the mass. Several dozen 
eclipsing binary systems have component mass and radius determinations that are good to 
1-2% (e.g., Andersen 1991), and show that main-sequence models for stars with masses 
in the range from about 1 Mq to 10 Mq and heavy element abundances near solar are in 
fairly good agreement with the observations. However, models for stars with masses that are 
significantly higher or lower, or that are in very early (pre-main sequence) or very advanced 
(post-main sequence) stages of evolution, or models for chemical compositions that are much 
different from solar, are largely untested by observations due to a lack of suitable systems or 
a lack of accuracy. 

Stellar evolution theory has a wide range of applications in modern astrophysics, some 
of which have profound cosmological implications. One such application is the estimate of 
the ages of globular clusters, which represent the oldest of the stellar populations in our 
Galaxy. Because the mctallicitics of globular clusters are typically much lower than solar, 
model fits to color-magnitude diagrams can be viewed as useful tests of theory at these 
metal-poor compositions, at least regarding the general shape of the isochrones. However, 
the predicted mass at any given point along a fitted isochrone cannot be tested directly 
against accurate observations because such a constraint is unavailable for low metallicities. 
While it is true that there is no immediate cause for concern regarding the accuracy of current 
metal-poor models, the use of such calculations for deriving the age and other properties of 
globular clusters remains an extrapolation to some degree, and it would be reassuring to 
have detailed observational support. 

Double- lined spectroscopic binary systems with known compositions below [Fc/H] = 
—0.5 that are also eclipsing and are therefore particularly suitable for accurate determinations 
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of the masses and radii are virtually non-existent, both in clusters and in the general field^. 
But the absolute masses can still be determined even when the binary is not eclipsing, if 
the pair can be spatially resolved so that the inclination angle of the orbit can be measured. 
With the highly precise astrometry now achievable using modern interferometers, this can 
significantly extend the pool of objects available for testing models of stellar evolution (e.g., 
Armstrong et al. 1992; Hummel et al. 1994; Boden et al. 1999a; Boden, Creech- Eakman, & 
Queloz 2000; Hummel et al. 2001). 

HD 195987 (HIP 101382, G209-35, Gliese 793.1, Groombridge 3215, a = 20^32™51':6, 

5 = +4r53'55", J2000. V = 7.08) is a binary with a period of 57.3 days (Imbcrt 1980) 
reported to have a mctallicity as low as [m/H]= —0.83 (Laird, Carney, & Latham 1988), 
although more recent estimates are not as low. Originally listed as a single-lined spectro- 
scopic binary, further observations of HD 195987 have revealed the weak secondary spectrum 
so that the system is now double-lined. The possibihty of eclipses has been mentioned oc- 
casionally in the literature, and would obviously make the system even more interesting. 
However, the astrometric observations described in this paper clearly rule that out. The 
object was marginally resolved on one occasion by the speckle technique (Blazit, Bonneau, 

6 Foy 1987), but subsequent attempts have failed to separate the components (Balega et 
al. 1999; Mason et al. 2001a). Measurements by the HIPPARCOS mission (Perryman et al. 
1997) detected the photocentric motion of the system with a semimajor axis of about 5 mas, 
and allowed the inclination angle and therefore the absolute masses to be measured for the 
first time, albeit with relatively large uncertainties (see Osborn & Hershey 1999). More 
recently the components have been clearly resolved by the Palomar Testbed Interferometer 
(PTI), opening the possibility of establishing the inclination angle of the orbit with much 
higher accuracy. This makes it a potentially important system in which the absolute masses 
can be determined and used to compare with stellar evolution theory for metal-poor stars, 
and is the motivation for this paper. 

We report here our spectroscopic and interferometric observations of HD 195987, the 
analysis of which has provided accurate masses (good to 1-2%) and luminosities for the 
two components. We compare these measurements with recent models for the appropriate 
met alii city. 



recent exception is the binary system OGLEGC-17 in the globular cluster iv Cen, discovered by the 
OGLE microlensing group (Udalski et al. 1993). A preliminary study by Thompson et al. (2001) yielded 
encouraging results with errors in the masses and radii of ~7% and ^3%, respectively, and a more detailed 
analysis using additional observations can be expected in the near future (sec Kaluzny et al. 2001). At an 
estimated metallicity of [Fe/H]= —2.29 this is an extremely important system. 
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2. Observations 
2.1. Spectroscopy 

HD 195987 was originally placed on the observing program at the Harvard-Smithsonian 
Center for Astrophysics (CfA) as part of a project to monitor the radial velocities of a sam- 
ple of high proper motion objects selected from the Lowell Proper Motion Survey (Giclas, 
Burnham, & Thomas 1971. 1978). The goal of this project is to investigate a variety of is- 
sues related to Galactic structure (see Carney & Latham 1987). Observations were obtained 
from 1983 July to 2001 January mostly with an echelle spectrograph mounted on the 1.5-m 
Wyeth reflector at the Oak Ridge Observatory (Harvard, Massachusetts), and occasionally 
also with nearly identical instruments on the 1.5-m Tillinghast reflector at the F. L. Whipple 
Observatory (Mt. Hopkins, Arizona) and the Multiple Mirror Telescope (also on Mt. Hop- 
kins, Arizona) prior to its conversion to a monolithic 6.5-m mirror. A single echelle order 
was recorded with photon-counting intensified Reticon detectors at a central wavelength of 
5187 A, with a spectral coverage of 45 A. The strongest features present in this window are 
the lines of the Mg I b triplet. The resolving power is A/AA ~ 35,000, and the signal-to-noise 
(S/N) ratios achieved range from about 13 to 50 per resolution element of 8.5 kms^^. A 
total of 73 usable spectra were obtained over a period of more than 17 years. The stability of 
the zero-point of our velocity system was monitored by means of exposures of the dusk and 
dawn sky, and small systematic run-to-run corrections were applied in the manner described 
by Latham (1992). 

The observations for the first 3 years were analyzed by Latham et al. (1988) with 
the techniques available at the time, but the lines of the secondary are so faint that they 
escaped detection and the object was treated as single-hned. With the introduction of more 
sophisticated analysis methods it was later discovered that the lines from the secondary 
could in fact be seen (Carney et al. 1994). A prehminary analysis has been presented by 
Goldberg et al. (2002) on the basis of the same spectra as the earher study. For the present 
investigation many more observations were obtained at higher S/N ratios and at phases that 
contribute more information to the mass determinations. Radial velocities from all of the 
spectra were derived by cross-correlation using the two-dimensional algorithm TODCOR 
(Zucker & Mazeh 1994). The templates used for the primary and secondary of HD 195987 
were taken from an extensive library of synthetic spectra based on model atmospheres by R. 
L. Kurucz^, computed by Jon Morse specifically for the wavelength window and resolution 
of our spectra. They are available for a range of effective temperatures (Teff), projected 
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rotational velocities (vsmi), surface gravities (log 51), and metallicities ([m/H]). Grids of 
correlations were run to establish the template parameters giving the best match to each 
component, based on the peak correlation value averaged over all exposures and weighted by 
the S/N ratio of each spectrum. The stars present no measurable rotational broadening in our 
spectra, so we adopted vsini values of kms~^, along with logg values of 4.5, appropriate 
for dwarfs. The temperature and metallicity determinations are complicated by the fact that 
the secondary component is extremely faint (only 10% of the light of the primary; see below), 
and in addition those two quantities are very strongly correlated. We fixed [m/H] at values 
ranging from —1.0 to +0.5 (in steps of 0.5 dex) and determined the effective temperatures 
in each case. Lower metaUicities lead to cooler temperatures because the changes in the 
strength of the lines from these two effects tend to compensate each other. Formally, the 
best match to the observed spectra was found for a metallicity around [m/H]= —0.2 and 
effective temperatures of = 5350 K and = 4550 K for the primary and secondary, 
respectively. However, because of the tradeoff mentioned above, different combinations of 
[m/H] and Tgs give nearly equally good fits without changing the velocities very much, and 
an external constraint is needed on one of these parameters to break the degeneracy. For this 
we chose the metallicity, since several determinations are available and appear to agree that 
the system is metal-deficient compared to the Sun (§5). We adopted the value [m/H]= —0.5, 
which leads to effective temperatures of = 5200 ± 150 K and = 4450 ± 250 K. The 
light ratio in our spectral window was found to be {Ib/Ia)5187 — 0.09 ± 0.01. 

Systematic errors in the radial velocities resulting from the narrow spectral window 
were investigated by means of numerical simulations, as discussed in detail by Latham et al. 
(1996). Briefly, we generated a set of artificial binary spectra by combining the primary and 
secondary templates in the appropriate ratio and applying velocity shifts for both components 
as computed from a preliminary orbital solution at the actual times of observation of each 
of our spectra. These artificial spectra were then processed with TODCOR in exactly the 
same way as the real observations, and the resulting velocities were compared with the 
input (synthetic) values. The differences were found to be smaller than 0.4 kms~^, but were 
nevertheless applied to the raw velocities as corrections. The effect on the minimum masses 
derived from the spectroscopic orbit, which depend on the velocity semi-amplitudes Ka and 
Kb, is at the level of 0.5%. The final velocities, including corrections, are listed in Table 1. 

The spectroscopic orbital elements resulting from these velocities are given in column 2 
of Table 2. They are the period (P, in days), the center-of-mass velocity (7, in kms~^), the 
radial velocity semi-amplitudes of the primary and secondary {Ka and Kb, in kms~^), the 
orbital eccentricity (e), the longitude of periastron of the primary {ua, degrees), and the time 
of periastron passage (T, Julian days). Earlier single-lined orbital solutions for HD 195987 
were published by Imbert (1980), Latham et al. (1988), and Duquennoy & Mayor (1991), 
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and Goldberg et al. (2002) recently reported a double-lined solution. They are included in 
Table 2 for comparison with ours. The orbit published by Latham et al. (1988), which is 
based on a small subset of the same spectra used in the present paper, was superseded by 
that by Goldberg et al. (2002) that used the same material, which in turn is superseded by 
our new definitive results. The solution by Duquennoy & Mayor (1991) is not independent 
of that by Imbert (1980), but represents an update using additional observations with the 
same instrument. As seen in the table, all these solutions are fairly similar. 

The uncertainties given in Table 2 for the velocity amplitudes of our new solution, upon 
which the masses of the components depend critically, arc strictly internal errors. In addition 
to the biases described earlier, further systematic errors in the velocities can occur because of 
uncertainties in the template parameters, particularly the temperatures and the metallicity. 
Extensive grids of correlations showed that the sensitivity of the velocity amplitudes to the 
secondary temperature is minimal because that star is so faint. They also showed that the 
combined effects of errors in the primary temperature and the metallicity contribute an 
additional 0.04 kms~^ uncertainty to Ka and 0.3 kms^^ to Kb- These have been combined 
with the internal errors in the amplitudes, and propagated through to the final masses that 
we report later in §5. 



2.2. Interferometry 

Near-infrared, long-baseline interferometric measurements of HD 195987 were conducted 
with the Palomar Testbed Interferometer (PTI), which is a 110-m baseline H- and fC-band 
(A ~1.6 ^m. and ~2.2 yum) interferometer located at Palomar Observatory. It is described in 
full detail elsewhere (Colavita et al. 1999). The instrument gives a minimum fringe spacing 
of about 4 mas at the sky position of our target, making the HD 195987 binary system 
readily resolvable. 

The interferometric observable used for these measurements is the normalized fringe 
contrast or visibility (squared), V^, of an observed brightness distribution on the sky. The 
analysis of such data in the context of a binary system is discussed in detail by Hummel et 
al. (1998), Boden et al. (1999a), Boden et al. (1999b), Boden, Greech-Eakman, & Queloz 
(2000), and Hummel et al. (2001), and will not be repeated here. 

HD 195987 was observed with PTI in conjunction with objects in our calibrator list in the 
K band (A ~2.2 //m) on 32 nights between 1999 June 24 and 2001 September 28, covering 
roughly 14 periods of the system. Additionally, the target was observed in the H band 
(A ~1.6 /xm) on five nights between 2000 July 20 and 2001 September 27. HD 195987, along 
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with the calibration objects, was observed multiple times during each of these nights, and 
each observation, or scan, was approximately 130 seconds long. For each scan we computed 
a mean value from the scan data, and the error in the V"^ estimate from the rms internal 
scatter (Colavita 1999). HD 195987 was always observed in combination with one or more 
calibration sources within ~20° on the sky. For our study we used four stars as calibrators: 
HD 195194, HD 200031, HD 177196, and HD 185395. Table 3 hsts the relevant physical 
parameters for these objects. The last two calibration objects are known to have at least 
one visual companion (Washington Double Star Catalog; Mason et al. 2001b) that could 
conceivable affect our visibility measurements. The companion of HD 177196 is so distant 
(44" at last measurement in 1925) that it has no effect because of the effective 1" field stop of 
the PTI fringe camera focal plane (Colavita et al. 1999). Of the three recorded companions 
of HD 185395 two are more than 40" away, and the other has an angular separation of 2"9. 
This close companion is more than 5 magnitudes fainter than the primary in the K band, 
and can therefore also be safely ignored. 

The calibration of the HD 195987 V"^ data was performed by estimating the interfer- 
ometer system visibility (Kjys) using calibration sources with model angular diameters, and 
then normalizing the raw HD 195987 visibility by V^^ to estimate the V"^ measured by an 
ideal interferometer at that epoch (Mozurkewich et al. 1991; Boden et al. 1998). Calibrating 
our HD 195987 dataset with respect to the four calibration objects listed in Table 3 results 
in a total of 171 calibrated scans (134 in K, 37 in H) on HD 195987 on 37 different nights. 
Our calibrated synthetic wide-band visibility measurements in the H and K bands are sum- 
marized in Table 4 and Table 5, which include the time of observation, the cahbrated V'^ 
measurement and its associated error, the residual from the final fit (see below), the {u,v) 
coordinates in units of the wavelength A (weighted by the S/N ratio), and the orbital phase 
for each of our observations. 



3. Determination of the orbit 

The radial velocities and interferometric visibilities for HD 195987 contain complemen- 
tary information on the orbit of the system, which is described by the 7 elements mentioned 
in §2.1 and the additional elements a (relative semimajor axis, expressed here in mas), i (in- 
clination angle), and Q (position angle of the ascending node). Only by combining data from 
both techniques can all 10 elements of the three-dimensional orbit be determined. While the 
center-of-mass velocity and the velocity amplitudes depend only on the spectroscopy, the 
information on a, i, and Q is contained solely in the interferometric visibilities. However, the 
V'^ values are invariant to a 180° change in the position angle of the binary, so distinguishing 
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between the descending node and the ascending node (where, according to convention, the 
secondary is receding) rehes on our knowledge of the radial velocities. 

As in previous analyses using PTI data (Boden et al. 1999a,b; Boden, Creech- Eakman, 
& Queloz 2000; Boden & Lane 2001), the orbit of HD 195987 is solved by fitting a Keplerian 
orbit model directly to the calibrated V'^ and radial velocity data simultaneously. This 
allows the most complete and efficient use of the observations''', so long as the two types of 
measurements are free from systematic errors (see below). It also takes advantage of the 
redundancy between spectroscopy and interferometry for the four elements P, e, uja, and T. 
Relative weights for the velocities (which are different for the primary and the secondary) 
and the visibilities were applied based on the internal errors of each type of observation. 

Formally the interferometric visibility observables have the potential for resolving not 
only the HD 195987 relative orbit, but the binary components themselves, and these two 
effects must be considered simultaneously (e.g.. Hummel et al. 1994). However, in the case of 
HD 195987 at a distance of approximately 22 pc the dwarf components have typical apparent 
sizes of less than 0.5 mas, which are highly unresolved by the 3-4 mas PTI fringe spacing. 
We have therefore estimated the apparent sizes of the components using the bolometric flux 
and effective temperatures (see Blackwell & Lynas-Gray 1994, and references therein), and 
constrained the orbital solutions to these model values (§5). We have adopted apparent 
component diameters of 0.419 ± 0.018 mas and 0.314 ± 0.049 mas for the primary and 
secondary components, respectively. These values are much smaller than the PTI H and K 
fringe spacings and have a negligible effect on the parameters of the orbit. 

The results of this joint fit are listed in Table 6 ("Pull-fit" solution), where we have 
used the //-band and T^'-band visibilities simultaneously, weighted appropriately by the 
corresponding errors. In addition to the orbital elements, the visibility measurements provide 
the intensity ratio in H and K, which are critical for the model comparisons described later. 

The uncertainties listed for the elements in Table 6 include both a statistical compo- 
nent (measurement error) and our best estimate of the contribution from systematic errors, 
of which the main components are: (1) uncertainties in the calibrator angular diameters 
(Table 3); (2) the uncertainty in the center-band operating wavelength (Aq ~ 1.6 /im and 



^In principle it would be possible to derive the position angle {9) and the angular separation (p) of the 
binary from the visibility measurements for each night. These intermediate data could then be used instead of 
the visibilities in what might be considered a more conventional astrometric-spectroscopic solution. However, 
the result would be inferior to the fit we describe here for several reasons (see, e.g., Hummel et al. 2001), not 
the least of which is the fact that, due to the pattern of our observations, the number of scans on many nights 
is insufficient to solve for the separation vectors, which would result in a rather serious loss of information 
(phase coverage). Therefore, we refrain from even listing values of {p, 6} in Table 4 and Table 5. 
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2.2 fj,m), taken to be 20 nm (~1%); (3) the geometrical uncertainty in our interferometric 
baseline (< 0.01%); and (4) listed uncertainties in the angular diameters assumed for the 
stars in HD 195987, that were held constant in the fitting procedure. 

Figure 1 depicts the relative visual orbit of the HD 195987 system on the plane of the 
sky, with the primary component rendered at the origin and the secondary component shown 
at periastron. We have indicated the phase coverage of our V"^ data on the relative orbit 
with heavy line segments. Our interferometric data sample essentially all phases of the orbit 
(see also Figure 4), leading to a reliable determination of the elements. The fit to these 
data is illustrated in Figure 2, which shows four consecutive nights of PTI visibility data on 
HD 195987 (24-27 June 1999), and predictions based on the "Full Fit" model for the 
system (Table 6). The fit to the radial velocity measurements is shown in Figure 3. Phase 
plots of the V"^ and velocity residuals are shown in Figure 4, along with the corresponding 
histograms. We note here in passing that the agreement between the speckle measurement by 
Blazit, Bonneau, & Foy (1987) (which by chance was obtained very near periastron passage, 
at phase 0.981) and our intcrferometric-spectroscopic orbit is very poor. However, this is 
not surprising given the limited resolution of the speckle observation (32 mas in the visible) 
compared to the predicted separation at the time of the observation (10.5 mas). 

For comparison, wc list also in Table 6 the solutions we obtain for HD 195987 using 
only the interferometric visibilities ("y^-only"), and only the radial velocities ("RV-only", 
repeated from Table 2). These two separate fits give rather similar results, indicating no 
significant systematic differences between the spectroscopic and interferometric data sets. 
Also included in the table is the solution reported in the HIPPARCOS Catalogue, where 
the semimajor axis refers to the photocentric motion of the pair rather than the relative 
motion. In the HIPPARCOS solution the elements P, e, loa, and T were adopted from the 
work by Duquennoy & Mayor (1991), and held fixed. Given the much larger formal errors, 
the resulting elements are as consistent with our solution as can be expected. 

4. The light ratio in the optical 

The interferometric measurements of HD 195987 with PTI provide the intensity ratio 
between the components in the H and K bands, and allow the individual luminosities to be 
determined in the infrared. A similar determination in the optical based on our spectroscopy 
was given in §2.1. A small correction from the 5187 A region to the visual band yields 
{Ib/La)v — 0.10 ± 0.02. Because the secondary is so faint, this estimate may be subject to 
systematic errors that are difficult to quantify. 
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The detection of the photocentric motion of the binary by HIPPARCOS, along with 
the measurement of the apparent separation with PTI, allows an independent estimate of 
the light ratio to be made in the optical. The relative semimajor axis (a) and the pho- 
tocentric semimajor axis (a) are related by the classical expression a = a{B — f3), where 
B = Mb/ {Ma + Mb) is the mass fraction (also expressed as -B = Ka/ {Ka + Kb) in terms 
of the observables) and (3 = Ib/{Ia + Ib) is the fractional luminosity. The ratio Ib/Ia at 
the effective wavelength of the HIPPARCOS observations {Hp passband) can therefore be 
determined. 

In order to take advantage of the much improved orbital elements compared to those 
available to the HIPPARCOS team, we have re-analyzed the HIPPARCOS intermediate 
astrometric data (abscissae residuals) as described in the original Catalogue (ESA 1997) 
(see also Pourbaix & Jorissen 2000). The orbital period, the eccentricity, the longitude of 
periastron, and the time of periastron passage were adopted from our combined visibility- 
radial velocity fit as listed in Table 6, and held fixed. In addition, the inclination angle and 
the position angle of the node are now known much better than can be determined from the 
HIPPARCOS data, so we adopted our own values here as well. The only parameters left to 
determine are then the semimajor axis of the photocenter (a), the corrections to the position 
and proper motion of the photocenter, and the correction to the HIPPARCOS parallax. 

The relevant results from this analysis are a = 5.04 ± 0.47 mas and tthip = 45.30 ± 
0.46 mas. The uncertainties may be somewhat underestimated because they do not account 
for correlations between the elements (since we have fixed several of them in the new solu- 
tion), although we do not expect the effect to be large in this case. As it turns out, the new 
results are not very different from the original determinations by the HIPPARCOS team 
{a = 5.24 ± 0.66 mas and ttrip = 44.99 ± 0.64 mas). Forcing the HIPPARCOS parallax to 
come out identical to our orbital parallax makes a negligible change in the semimajor axis 
of the photocenter. 

The luminosity ratio we determine from this is {Ib/Ia)v ~ 0.13± 0.04, which includes a 
very small correction from the band to the Johnson V band. This is in good agreement 
with our spectroscopic determination, and we adopt the weighted average of the two results, 
{Ib/Ia)v = 0.11 ± 0.02, which corresponds to = 2.4 ± 0.2 mag. 

5. Physical properties of HD 195987 

The absolute masses of the components follow from the results of our combined interf- 
erometric-spectroscopic solution summarized in Table 6, and are Ma — 0.844 ± 0.018 Mq 
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and Mb — 0.6650±0.0079 Mq^. The high accuracy achieved (2.1% error for the primary and 
1.2% error for the secondary) is the result of the high quahty of the observations as well as 
the favorable geometry of the system (sini ^ 0.99). The limiting factor is the spectroscopy, 
and in particular the velocities for the faint secondary. 

The combination of the angular separation from interferometry and the linear separation 
from spectroscopy yields the orbital parallax of the system independently of any assumptions 
beyond Newtonian physics. The result, TTorb = 46.08 ± 0.27 mas (distance= 21.70 ± 0.13 pc) 
is nearly a factor of two more precise than our revised HIPPARCOS value of tthip = 45.30 ± 
0.46 mas, and the two determinations differ by about 1.5a (1.7%). The corresponding 
difference in the distance modulus, and consequently in the absolute magnitudes derived 
below, is A(m — M) — 0.037 mag. The ground-based trigonometric parallax of HD 195987 
from the weighted average of 5 determinations is tt = 49.1 ± 5.1 mas (van Altena, Lee, & 
Hoffleit 1995). 

The component luminosities in the V, H, and K passbands depend on the flux ratios 
that we have measured, on the orbital parallax, and on the combined-light photometry. 
For the V magnitude of the system we adopt the average value given by Mermilliod & 
MermiUiod (1994), which is the result of 10 individual measurements from 7 different sources: 
V = 7.080 ± 0.016. Two measurements of K are available from Voelcker (1975) and from 
the 2MASS Catalog (Skrutskic ct al. 1997). Conversion of both to the CIT system (Elias 
et al. 1983), using the transformations by Besscll & Brett (1988) and Carpenter (2001), 
gives K = 4.98 ± 0.03. The H magnitude of the system measured by Voelcker (1975) 
is H = 5.23 ± 0.05. A 2MASS measurement in H is not available because the object is 
saturated in this passband, although not in J and K. This is somewhat unexpected if we 
assume normal colors for stars of this mass, and it casts some doubt on the reliability of 
the H-hand measurement. As a check we obtained a new measurement using an infrared 
camera on the 1.2-m telescope at the F. L. Whipple Observatory equipped with a 256 x 256 
InSb detector array, with standard stars adopted from Elias et al. (1982) and 2MASS. The 
result, H = 5.04 ± 0.03 (on the CIT system), is significantly different from the previous 
estimate, confirming our suspicions, and it is the value we adopt for the remainder of the 
paper. Any extinction and reddening corrections are negligible at a distance of only 22 pc. 
The individual absolute magnitudes and colors are listed in Table 7 along with the other 
physical properties. 

The V — K color indices for each component allow for an independent estimate of the 
effective temperatures of the stars. Based on the color-temperature calibrations by Martmez- 



^These uncertainties include the contribution due to the template parameters discussed in §2.1. 



- 12 - 



Roger, Arribas, & Alonso (1992), Alonso, Arribas, & Martmez-Roger (1996), and Carney et 
al. (1994) (the later being for the CIT system), we obtain average values of = 5200±100 K 
and — 4100 ± 200 K that are nearly independent of the adopted metal abundance (for 
which these calibrations include a corrective term). The primary estimate is identical to 
our spectroscopic value in §2.1, while the secondary estimate is somewhat lower than the 
spectroscopic value and its uncertainty has a substantial contribution from the photometric 
errors. 

An additional estimate of the effective temperature may be obtained from photometry 
in other bands, via dcconvolution of the combined light. For this we used the indices for 
HD 195987 in the Johnson and Stromgren systems as listed by Mermilliod & Mermilliod 
(1994) {B-V = 0.800 ± 0.014; U-B = 0.370 ± 0.026) and by Olsen (1983, 1993) {b-y = 
0.480 ± 0.002, mi = 0.296 ± 0.005, ci = 0.271 ± 0.005). The deconvolution was performed 
using tables of standard colors for normal stars by Lejeune, Cuisinier, & Buser (1998) (for 
a metallicity of [Fe/H]= —0.5) and by Olsen (1984), and adopting a magnitude difference 
between the primary and secondary in the visual band of AV = 2.4 mag (§4). In addition 
to the calibrations mentioned above we used those by Carney (1983) and Olsen (1984). 
The results based on the deconvolved B — V, b — y, and also the V — K indices for the 
system (accounting for the small difference between the different photometric systems in the 
infrared) give values of 5200 ± 100 K for the primary and 4100 ± 200 K for the secondary, 
identical to the estimates from the observed V — K colors. 

We adopt for the stars in HD 195987 the weighted average of our three determinations 
of the effective temperature (one spectroscopic and two photometric): Tjp = 5200 ± 100 K 
and Tgg = 4200 ±200 K. Prom these effective temperatures, an estimated system bolometric 
flux of (5.308 ± 0.062) x 10~^ erg cm~^ s~^ based on archival broadband photometry, and 
our observed 2.2 iim intensity ratio we estimate component bolometric fluxes of (4.28±0.19) 
and (1.02 ± 0.28) x 10~^ erg cm~^ s~^. These combined again with the adopted component 
temperatures yield estimated apparent diameters of 0.419 ±0.018 mas and 0.314 ±0.049 mas 
for the primary and secondary components, respectively, which are the values adopted for 
the orbital fit (§3). At our estimated system distance to HD 195987 these angular diameters 
imply physical radii of 0.979 ± 0.043 R© and 0.73 ± 0.11 R© for the primary and secondary 
components, respectively. 

As mentioned in §1, early high- resolution spectroscopic estimates of the metal abun- 
dance of HD 195987 placed the system in an interesting range for absolute mass determi- 
nations, at a metallicity of [m/H]= —0.83 ± 0.15 (Laird, Carney, & Latham 1988). Subse- 
quently the same authors revised their estimate upward to [m/H]= —0.60 ± 0.15 using the 
same techniques and more spectra (Carney et al. 1994). Other determinations in the litera- 
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ture have given values that are similar or somewhat closer to the solar abundance. Fulbright 
(2000) used high-resolution spectroscopy to derive an estimate of [Fe/H]= —0.66 ± 0.13, 
while Beers et al. (1999) obtained [Fe/H]= —0.52 ± 0.20 with a combination of different 
spectroscopic methods based on lower-resolution spectra. Marsakov & Shevelev (1988) in- 
ferred [Fe/H]= —0.34 based on the ultraviolet excess of the system, and Wyse & Gilmore 
(1995) used Stromgren photometry to estimate the metallicity at [Fe/H]= —0.31. Despite 
the range of values, these sources do seem to point toward a heavy element abundance mod- 
erately lower than the Sun, by perhaps a factor of 3 or so. We discuss possible biases in the 
metallicity estimates in §6.2. 

6. Discussion 

With many of the fundamental physical properties of HD 195987 now known to high 
accuracy, we proceed in this section with a detailed comparison between the observations and 
recent models of stellar evolution. In addition to the prospect of an interesting comparison 
for a chemical composition lower than solar, both components of our system have masses 
below 1 Mq, a regime in which relatively few stars have absolute mass determinations good 
to 1-2% (see, e.g., Andersen 1991; Clausen et al. 1999; Delfosse et al. 2000). 

6.1. Comparison with stellar evolution models 

In order for a given model to be successful, we require that it agree with all the mea- 
surements for HD 195987 simultaneously, and for both components at the same time since 
they are presumably coeval. The measurements are the absolute masses, the absolute mag- 
nitudes in three different passbands (My, Mh, Mk), the effective temperatures, and the 
metal abundance. 

In Figure 5 we show the absolute magnitudes and temperatures of the primary and 

secondary of HD 195987 as a function of mass, compared with isochrones from the Yale- 
Yonsei models by Yi et al. (2001)^. Changes in the isochrones with age for a fixed metallicity 
equal to solar are illustrated on the left, and consist not only of a vertical displacement but 
also of a change in slope, becoming steeper for older ages. A similar effect is seen in the 



^The infrared magnitudes in these isochrones are based on the color tables by Lejeune, Cuisinier, & 
Buser (1998), which adopt filter transmission functions from Bessell & Brett (1988) for the Johnson-Glass 
system. For the comparison with our observations we transformed the isochrones to the CIT system using 
the corrections by Bessell & Brett (1988). 
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diagrams on the right-hand side, which show changes in the models for a fixed age (12 Gyr) 
and a range of metalhcities bracketing the values reported for HD 195987. The isochrones 
shown with a solid line, which correspond to solar metallicity and an age of 12 Gyr, provide 
a fairly good fit to the absolute magnitudes in the H and K bands, as well as to the My of 
the primary star. However, the secondary component appears undcrluminous in V , or else 
the slope of the models in the mass- My plane is too shallow. In addition, the temperatures 
predicted by the models may be somewhat too hot, particularly for the secondary. Setting 
these differences aside for the moment, the Yi et al. (2001) models would appear to point 
toward a metallicity for HD 195987 that is not far from solar, along with a fairly old age. This 
metallicity seems to confiict with the observational evidence indicated earlier. Formally the 
best fit to the magnitudes in H, K, and V (primary only) gives a metallicity of [Fe/H]= —0.07 
and an age of 11.5 Gyr. Lowering the age of the isochrones to 10 Gyr and at the same 
time decreasing the metallicity to [Fe/H]= —0.15 still produces tolerably good fits, but 
metalhcities much lower than this are inconsistent with the measured magnitudes for any 
age. The disagreement with the My of the secondary and with the temperatures remains. 

In order to investigate the discrepancy in the absolute visual magnitude for the sec- 
ondary, we focus in Figure 6 on the mass- My plane and add observations for other binary 
components that have accurately determined masses and luminosities: FL Lyr B, HS Aur A, 
and HS Aur B (Andersen 1991); Gliese 570 B and Gliese 702 B (Delfosse et al. 2000); and 
V818 Tau B and YY Gem AB (Torres & Ribas 2002). Figure 6a shows that it is not only the 
secondary of HD 195987 (filled circles) that appears undcrluminous compared to the Yi et al. 
(2001) isochrones, but four other stars with similar masses suggest the same trend (triangles). 
In Figure 6b we show the same observations compared to the Lyon models by Baraffe et al. 
(1998), for the same ages and heavy-element composition as in Figure 6a. The agreement 
for the secondary is now much better, and is presumably due to the use by Baraffe et al. 
(1998) of sophisticated model atmospheres as boundary conditions to the interior equations, 
whereas the Yi et al. (2001) models use a gray approximation. The latter has been shown to 
be inadequate for low mass stars, where molecular opacity becomes important (Chabrier & 
Baraffe 1997). Figure 6 suggests that the discrepancy begins at a mass intermediate between 
that of the primary and the secondary in HD 195987. One other consequence of the use of 
the gray approximation is that the effective temperature is typically overestimated in those 
models, just as hinted by Figure 5. 

In Figure 7 we illustrate the agreement between the observations for HD 195987 and the 
predictions from the models by Baraffe et al. (1998) and also by Siess, Forestini, & Dougados 
(1997) (see also Siess, Dufour, & Forestini 2000). The latter isochrones also use a non-gray 
approximation (from fits to a different set of model atmospheres than those used by the Lyon 
group), and should therefore be fairly realistic as well. Once again the Baraffe et al. (1998) 
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models (left-hand side) show reasonably good agreement in the mass-luminosity planes for 
solar composition and an old age (12 Gyr), but do not quite reproduce the estimated effective 
temperature of the primary. The Siess, Forestini, & Dougados (1997) models are similar, 
although the agreement in My is not as good and the fit to the primary temperature is 
somewhat better. 

In addition to the boundary conditions, another difference between the Baraffc ct al. 
(1998) models and the Yi ct al. (2001) models is the extent of mixing allowed, as described 
in the standard convection prescription by the mixing-length parameter Oml- The models 
by Yi et al. (2001) adopt a value of q;ml = 1-74 that best fits the observed properties of the 
Sun. A much lower value of q;ml = 1-0 is used by Baraffe et al. (1998), whereas the best 
fit to the Sun in those models requires Q!ml = 1-9. This has significant consequences for the 
temperature profile and other properties such as the radius. In particular, less mixing leads 
to a lower effective temperature. The effect is illustrated in the diagrams on the left side 
of Figure 7, where the dot-dashed lines gives the predictions from the Baraffe et al. (1998) 
models for [Fe/H]= 0.0 and an age of 12 Gyr if the solar value of ckml is used. Changes in 
Mki Mh, and My compared to the isochrones for the same age and metallicity but with 
c^ML = 1-0 (solid lines) are relatively small and in fact tend to improve the agreement in My, 
while the temperatures predicted with qjml = 1-9 are several hundred degrees hotter than 
with ctML = 1.0 and show better agreement for the primary (the slope is also a better match 
to the observations). Presumably an intermediate value for the mixing- length parameter 
used in the Baraffe et al. (1998) models would provide an optimal fit to HD 195987 for the 
composition and age indicated. 

Other recent models give similar fits to the observations for HD 195987. For example, 
the recent series of isochrones from the Padova group (Girardi et al. 2000) give fits that are 
also good for a metaUicity near solar, and a slightly older age of 14 Gyr. 

A shortcoming of the comparisons above is that a number of additional parameters in 
the isochrones are fixed, and they are somewhat different for each series of models. One 
of such parameters is the helium abundance, Y. The calculations by Yi et al. (2001) use 
Y = 0.266 for solar metallicity (Figure 5), while Baraffe et al. (1998) adopt Y = 0.275 for 
their models with auL = 1-0, but use a higher value of F = 0.282 for the more realistic 
models that fit the Sun with auL = 1-9 (Figure 7). Siess, Forestini, & Dougados (1997) 
adopted Y = 0.277, and Girardi et al. (2000) used Y = 0.273. While the differences in these 
adopted helium abundances for solar metallicity are not large, they do affect the comparisons 
to some degree. An increase in Y will shift the isochrones upwards in the mass-luminosity 
diagrams (opposite effect as an increase in metallicity, Z), and will yield slightly higher 
effective temperatures. Changes in Y for other metallicities are governed by the enrichment 
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law adopted in each series of models. 

The effect of the treatment of convective overshooting is illustrated on the right-hand 
side of Figure 7 for the Siess, Forestini, & Dougados (1997) models. The dot-dashed hne 
corresponds to the calculations for an age of 12 Gyr and an overshooting parameter of 
CKov = 0.2Hp (where Hp is the pressure scale height), while the other isochrones assume no 
overshooting. The change compared to the 12 Gyr isochrone with no overshooting (solid 
line) is hardly noticeable. 

6.2. The iron abundance of HD 195987 

The model comparisons above seem to point toward a metallicity for the system near 
solar that is at odds with the observations, along with a fairly old age (~10-12 Gyr), perhaps 
a somewhat unusual combination. Because of the possible implications of this disagreement 
for our confidence in the models, in this section we examine each metallicity determination 
more closely in an effort to understand these discrepancies, and we attempt to quantify 
possible systematic errors. 

The presence of the secondary, even though it is faint, affects the total light of the 
system at some level introducing subtle biases in the photometric estimates of the metallicity 
of HD 195987, which are derived from the combined colors assuming that they correspond to 
a single star. The photometric estimates of the temperature based on the same assumption 
are also affected, and this may propagate through and affect some of the spectroscopic 
abundance determinations as well. 

In order to quantify these effects we have simulated binary systems by combining the 
single-star photometry of a primary and a secondary, each with normal colors, based on the 
same tabulations used in §5. We computed various photometric indices from the combined 
hght, and used them to estimate the metaUicity of HD 195987 following the same procedures 
employed by Marsakov & Shevelev (1988) (who relied on the ultraviolet excess in the Johnson 
system) and Wyse & Gilmore (1995) (based on Stromgren photometry). We then compared 
these results with those obtained for the primary alone, for a range of magnitude differences 
AV between the primary and secondary. In all cases we selected the primary so as to 
reproduce the actual observed colors of HD 195987 in each photometric system at AV = 
2.4 mag. 

Figure 8a shows how the presence of the secondary affects the U — B and B — V colors, 
as well the ultraviolet excess S{U — B) used by Marsakov & Shevelev (1988), as a function 
of the magnitude difference. The normalized ultraviolet excess 6{U—B)q q (corrected for the 
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guillotine; see Sandage 1969), also frequently used for metallicity estimates, has a similar 
behavior. Figure 8b illustrates the effect on the derived metallicity using the calibration by 
Marsakov & Shevelev (1988), and also two calibrations by Carney (1979) based on d{U—B)Q Q. 
The [Fe/H] estimates are biased toward lower values^'^, and the maximum effect is for a 
secondary about 2 magnitudes fainter than the primary. If the original determination by 
Marsakov & Shevelev (1988) ([Fe/H]= -0.34) is corrected for this effect (-0.2 dex at AV = 
2.4 mag), the result is much closer to the solar abundance. 

The normalized ultraviolet excess observed for HD 195987 is 6{U — B)o.e = 0.061. This 
leads to [Fe/H]= —0.13 or [Fe/H]= —0.18 using the linear or quadratic calibration formulae 
by Carney (1979). From Figure 8b it is seen that in this case the corrected values would be 
[Fe/H]~ +0.1. 

Similar effects occur in the Stromgren system. The changes in the photometric indices 
are shown in Figure 9a. The estimate of [Fe/H]= —0.31 by Wyse & Gilmore (1995) is 
based on the quantities {b — y, mi, ci} and the calibration by Schuster & Nissen (1989). A 
different calibration by Olsen (1984) involving {5mi, Sci} yields a nearly identical value of 
[Fe/H]= —0.32. As seen in Figure 9b the bias toward lower metaUicities due to the secondary 
happens to reach a maximum quite near the AV of our system. If accounted for, the resulting 
metal abundance from the Stromgren indices would be essentially solar. 

Though it may seem that these corrections bring the photometric estimates of the 
metaUicity of HD 195987 in hne with the indications from the model comparisons described 
earher, the photometric determinations of [Fe/H] are as a rule less trustworthy than those 
based on high- resolution spectroscopy (see, e.g., Gehren 1988). 

The spectroscopic determinations in turn are not, however, completely insensitive to 
the presence of the secondary either. Two different effects must be considered. On the one 
hand the continuum from the secondary typically tends to fill in the spectral lines of the 
primary, which then appear weaker as if the star were more metal-poor. On the other hand 
the combined-light photometry is reddened. Therefore, the temperature estimates that are 
used to begin the spectroscopic analysis are biased toward lower values. 



^°The sign of these changes may seem counterintuitive. The secondary star tends to redden both the 
U—B and the B—V indices, which ought to make the combined Hght of the system appear more metal-rich. 
As it turns out, the slope of this "reddening" vector in the U — B vs. B — V plane is smaller than that of 
the blanketing vector, so that the net effect is a bias toward larger ultraviolet excesses, and therefore lower 
metallicity estimates. We find, however, that in general this depends both on the color of the primary and 
on the magnitude difference AV. For bluer primaries than that of HD 195987 the sign of the change is 
reversed in certain ranges of AV. 
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The metallicity determination by Carney et al. (1994), [m/H]= —0.60, which super- 
sedes that by Laird, Carney, & Latham (1988), is based on a subset of the same spectra 
that we have used in this paper, which are of high resolution but low S/N. To derive the 
metal abundance they used a technique to compare the observations against a grid of 
synthetic spectra for a range of metallicities, adopting a fixed temperature determined from 
photometric indices. As they discuss, the [m/H] values derived for double-lined binaries 
with this method can be quite different depending on whether the spectra that are used 
have the lines of the two components well separated or exactly aligned. In the first case the 
continuum from the secondary tends to fill in the lines of the primary, as mentioned earlier, 
and the metallicity is biased toward lower values. In the second case, when the velocities of 
the stars are similar (close to the center-of-mass velocity 7 of the binary), the effect is the 
opposite because the average line strength is slightly increased due to the contribution of the 
cooler secondary. Carney et al. (1994) used only the few spectra near the 7 velocity available 
to them, because then the bias toward higher metallicities is balanced to some degree by 
the bias toward lower metallicities resulting from the lower temperature inferred from the 
combined fight. The latter is the result of the strong correlation between metallicity and 
temperature, which is essentially the same as that mentioned in §2.1. 

In Figure 10a we illustrate the magnitude of the effect of the secondary on the broad- 
band and intermediate-band color indices B — V, b — y, and also V — K, from simulations 
analogous to those performed earlier. The contamination is much more noticeable in V—K, 
of course, because the secondary becomes comparatively brighter in K. Changes in the 
effective temperature inferred from the combined light compared to those for the primary 
alone are shown in Figure 10b, based on the same calibrations used by Carney et al. (1994). 

Simulations performed by Carney et al. (1994) to estimate the residual bias on the 
metallicity for double-lined systems indicate that in general the two effects do not quite 
cancel each other out. We estimate from their experiments that their metal abundance 
[m/H] for HD 195987 may still be too low by 0.15-0.20 dex, accounting for the fact that this 
system is a bit more extreme than the simulated binaries they considered. Note also that this 
metallicity estimate is labeled [m/H] instead of [Fe/H] because it includes all metals with 
lines present in the spectral window. The distinction could in principle be an important 
one due to the presence of the strong Mg I b lines in the CfA spectra, and the typically 
enhanced abundance of this element relative to iron in metal-poor stars (see next section). 
Tests described by Carney et al. (1987) indicate that the correction required for this effect is 
not likely to be more than about —0.05 dex for this system, especially given that it may be 
partially masked by other corrections they applied to the abundances based on comparisons 
with metallicity standards that may also have enhanced magnesium (J. Laird & B. Carney 
2002, priv. comm.). The contribution of the two effects described above leads to an adjusted 
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Carney et al. (1994) estimate of [Fe/H] between -0.50 and -0.45. 

The metallicity determination by Beers et al. (1999), [Fe/H]= —0.52, is based on lower- 
resolution spectra (1-2 A) and relies on both the strength of a Ca II K-line index and on the 
height of the peak of the Fourier autocorrelation function of the spectrum, which is a line 
strength indicator. Both indices are calibrated against the observed B—V so as to remove the 
temperature dependence. The unrecognized presence of the secondary is unlikely to weaken 
the spectral lines of the primary very much in this case since the observations cover only the 
blue portion of the spectrum, where the secondary is fainter. However the reddening in the 
B — V color, estimated to be ~0.04 mag from Figure 8a, may have some effect on [Fc/H]. 
It is difficult to quantify its effect due to the complexity of the technique used by Beers et 
al. (1999), but on the basis of their Figure 6 and Figure 9 the bias will be toward lower 
abundances for both of their metallicity indices, and a correction of +0.1 dex or more may 
be appropriate. With this adjustment the value for HD 195987 would be [Fe/H]f=i —0.45 to 
-0.40. 

Finally, the detailed high-resolution (A/AA ~ 50,000) abundance analysis by Fulbright 
(2000) gives [Fe/H]= —0.66, but also does not account for the presence of the secondary 
star. The measured equivalent widths and consequently the metal abundances for iron and 

other elements are thus expected to be slightly underestimated due to light of the secondary 
filling in the lines of the primary, given that their spectrum was obtained at an orbital phase 
when the lines were not exactly aligned. The magnitude of the effect may depend on the 
strength of the lines used in the analysis, and could be as large as 0.1 dex. In addition, 
Fulbright (2000) adopted an initial estimate of the effective temperature based on the V—K 
index, which according to Figure 10b leads to an underestimate of ATes ~180 K due to 
the infrared excess. Although their use of iron lines of different strength (to constrain the 
microturbulence) , different excitation potentials (to constrain Teflf), and different ionization 
stages (to contrain log g) should tend to reduce the sensitivity of the results to this problem, 
the parameter sensitivity experiments in Table 8 by Fulbright (2000) indicate that such an 
error in Tcs can still produce changes in the resulting abundance of 0.07-0.16 dex (depending 
on whether other parameters are held fixed or allowed to vary to compensate for ATeg). The 
changes will be in the direction of lower abundances for lower temperatures. The combined 
effect of these corrections results in a metallicity of [Fe/H]Ri —0.50 to —0.40. 

It may seem that many of these adjustments to account for the secondary are rather 
small and perhaps even within the uncertainties of some of the determinations. Indeed 
they are, but they are also systematic in nature and known to be present. And because 
we are trying to understand a systematic difference between the estimates of [Fe/H] and 
indications from all the models, these biases cannot be ignored, particularly since they often 
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go in the same direction. To summarize the state of the (corrected) metaUicity estimates 
for HD 195987, it appears that both of the photometric determinations result in abundances 
very near solar, while the three spectroscopic estimates, which should be more reliable, yield 
[Fe/H] values between —0.50 and —0.40. These are still decidedly lower than the solar value, 
and so the apparent disagreement with the evolutionary models persists if we give preference 
to the spectroscopic abundances. 

6.3. Detailed abundances for other elements; further model comparisons 

The detailed chemical analysis by Fulbright (2000) has revealed certain patterns in 
HD 195987 that are quite typical in metal-poor stars, and that turn out to be key to un- 
derstanding the comparison with theoretical isochrones. In particular, it is well known that 
the abundance of the so-called a-elements (O, Ne, Mg, Si, S, Ar, Ca, and Ti) in metal-poor 
stars is usually enhanced relative to iron when compared to the same ratios in Sun (Conti 
et al. 1967; Greenstein 1970). The enhancement, [a/Fe], is actually observed to depend on 
the metaUicity as a result of enrichment from Type II and Type la supernovae on differ- 
ent timescales throughout the history of the Galaxy (see, e.g.. Wheeler, Sneden, & Truran 
1989). The average enhancements for the different a-elements increase from the solar val- 
ues at [Fe/H]= 0.0 to about +0.3 to +0.5 dex at [Fe/H]~ —1 or so, and then they remain 
approximately constant for more mctal-dcficient stars. There is even evidence that different 
stellar populations in the Galaxy are chemically distinct (see, e.g., Edvardsson et al. 1993; 
Mashonkina & Gehren 2000; Prochaska et al. 2000). Four of the a-elements in HD 195987 
have been measured by Fulbright (2000), and have been found to be enhanced. The ratios 
measured are [Mg/Fe]= +0.44 ± 0.07, [Si/Fe]= +0.41 ± 0.07, [Ca/Fe]= +0.26 ± 0.07, and 
[Ti/Fe]= +0.31 ± 0.09 in the standard logarithmic notation. 

A stellar mixture rich in these elements produces an increase in the opacity that affects 
the structure of a star and changes both the temperature and the luminosity predicted by the 
models. Compared to model calculations with solar ratios for all elements, ct-enhancement 
leads to cooler temperatures and lower luminosities. Therefore, these anomalies must be 
accounted for if a proper interpretation is to be made of the observations. A number of 
authors have noted, however, that these effects can be mimicked in models that assume 
solar ratios by simply increasing the overall metaUicity used to compute the tracks (see, 
e.g., Chicffi, Straniero, & Salaris 1991; Chaboycr, Sarajcdini, & Demarque 1992; Baraffe et 
al. 1997; VandenBerg et al. 2000). Prescriptions for how to do this as a function of the a- 
enhancement have been presented by Salaris, Chieffi, & Straniero (1993) and VandenBerg et 
al. (2000), among others. In the absence of a detailed measurement of the abundance of the a- 
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elements for a particular object, average values of [ct/Fe] have often been adopted in different 
ranges of [Fe/H]. This carries some risk, however, given that there is considerable scatter in 
the measured enhancements for different stars at any given [Fe/H] value (Edvardsson et al. 
1993; Carney 1996; Prochaska et al. 2000). 

In the case of HD 195987 we have a direct measurement of [a/Fe] if we assume that 
all the a-elements follow the trend of Mg, Si. Ca, and Ti. which were the ones actually 
measured. This is generally observed to be true to first order in metal-poor stars. The 
average overabundance for the four measured elements is [Q;/Fe]= +0.36 ± 0.12^^. Following 
VandcnBcrg et al. (2000), the overall mctallicity adjustment required for this level of [a/Fc] is 
+0.27 dex. Given the estimates of [Fc/H] in the range —0.50 to —0.40 (§6.2), one would then 
expect that model isochrones computed for a metallicity between —0.23 and —0.13 should 
provide reasonably good fits to the observations. This is in fact quite close to what we find, 
as shown by the comparisons in §6.1, which suggested an overall composition between solar 
and -0.15. 

Though adjusting the metallicity of models with scaled-solar mixtures may mimic those 
with a-enhancements to a good approximation, the match is not perfect, as shown by Van- 
denBerg et al. (2000). Isochrones incorporating a-enrichment for a range of values of [a/Fe] 
between 0.0 to +0.6 have recently been published by Bcrgbusch & VandcnBerg (2001). We 
compare them with the observations for HD 195987 in Figure 11, in the My vs. mass and 
Teff vs. mass diagrams {H- and i^-magnitudc predictions are not available in this series of 
calculations). Because these models use a gray approximation for the boundary conditions 
between the photosphere and the interior, we do not expect them to reproduce the brightness 
of the secondary or the temperature of either star very well (see §6.1). Thus we must rely 
only on the absolute visual magnitude of the primary component for this test. 

For a fixed value of [a/Fc]= +0.36 the iron abundance that provides the best match is 
[Fe/H]= —0.45 for an age of 12 Gyr (Figure 11), and [Fe/H]= —0.52 for 10 Gyr. These values 
are in excellent agreement with the measured abundances by Carney et al. (1994), Beers et al. 
(1999), and Fulbright (2000), after correcting for the biases due to the secondary as described 
before. The effect of the a-enhancement is illustrated in the figure by the shift between the 



Because we are dealing in this case with ratios of element abundances, the effect that the secondary 

component might have on those determinations is not expected to be as significant as in the case of [Fe/H]. 
Table 8 by Fulbright (2000) indicates that the bias due to temperature errors averaged over Mg, Si, Ca, and 
Ti (which have different signs for the correlation with Tes) is 0.03 dex at most, with the sign depending on 
whether only the temperature is changed or whether other fitted parameters are allowed to vary simulta- 
neously to compensate. Similarly, the continuum from the secondary should have a minimal effect on the 
abundance ratios. 
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solid line ([a/Fe]^ +0.36) and the dashed hne ([o;/Fe]= 0.0) for the same metaUicity. The 
uncertainty in the adopted value of [a/Fe] (0.12 dex) translates into the same uncertainty 
for the fitted value of [Fe/H]. The discrepancy for the My of the secondary and for the 
temperatures is in the same direction as shown earlier, i.e., the models arc ovcrluminous for 
a star with the mass of the secondary, and both temperatures are overestimated. 

In addition to the a-elements, abundance ratios relative to iron were measured by Ful- 
bright (2000) for the light elements Na and Al, the iron-peak elements V, Cr, and Ni, and 
the heavy elements Ba, Y, and Eu. A comparison with similar measurements for other 
stars reported by Prochaska et al. (2000) reveals that the pattern of enhancements seen in 
HD 195987 is very similar to that exhibited by thick disk stars, and is clearly distinct from 
that of the thin disk stars (see also Fuhrmann 1998). On that basis one might conclude 
that HD 195987 belongs to the thick disk population of our Galaxy. This would typically 
imply a rather old age for the binary, in qualitative agreement with the indications from the 
models in §6.1. The kinematics of the system computed from the proper motion, the sys- 
temic radial velocity, and our orbital parallax (Galactic velocity components [/ = 18 kms~"^, 
y = 13 kms^^, W = 4& kms~^, relative to the Local Standard of Rest) lend some support 
to this idea: the W velocity is quite consistent with thick disk membership (see Garney, 
Latham, & Laird 1989), although the V component is not as extreme. 

The abundances of Eu and Ba are of particular astrophysical interest. Europium in 
the solar system is produced mostly by rapid neutron capture ("r-process"), for which the 
principal formation site is believed to be Type 11 supernovae. It is one of the few r-process 
elements with relatively unblended atomic lines present in the visible part of the spectrum, 
which makes it very useful as an indicator for the r-process history of stellar material. 
Barium, on the other hand, is predominantly a slow neutron capture ( "s- process" ) element 
that is thought to be produced mainly during the thermal pulses of asymptotic giant branch 
stars. The abundance ratio between Eu and Ba is therefore particularly sensitive to whether 
nucleosynthesis of heavy elements occured by the s- or r-process, and as such it is a useful 
diagnostic for studying the chemical evolution of the Galaxy. [Eu/Ba] is observed to increase 
with metal deficiency, and has also been shown to correlate quite well with age (Woolf, 
Tomkin, & Lambert 1995). For HD 195987 Fulbright (2000) measured [Eu/Fe]= +0.29±0.10 
and [Ba/Fe]= -0.10 ±0.12, so that [Eu/Ba]= +0.39 ±0.14. Once again this is quite typical 
of the values for thick disk stars (see Mashonkina & Gehren 2000). Furthermore, according 
to the correlation reported in Figure 14 by Woolf, Tomkin, & Lambert (1995) this ratio 
suggests an age for HD 195987 that is at least 10-12 Gyr, in agreement with indications 
from the stellar evolution models used here and also consistent with the old age of the thick 
disk (Garney, Latham, & Laird 1989). 
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6.4. The overall agreement with the models 

The comparisons in §6.1 and §6.3 indicate that current stellar evolution theory can 
successfully predict all of the observed characteristics of HD 195987 to high precision for 
a chemical composition matching what is measured. Unfortunately, however, none of the 
published models is capable of giving a good fit to all properties of both stars simultaneously. 
The various isochrones used in this paper do at least agree in predicting an old age for the 
system, which is also consistent with the pattern of enhancement of the heavy elements, as 
described above. 

The low luminosity of the secondary star in the visual band lends strong support to 
models that use sophisticated boundary conditions to the interior equations, as in Baraffe 
et al. (1998), particularly for the lower main sequence. This also appears to be required in 
order to reproduce the effective temperatures of these stars, but not with a mixing-length 
parameter as low as in those models. A value of Oml much closer to the one for the Sun (as 
used in most models except those by Baraffe et al. 1998) seems indicated. 

The enhanced levels of the a-elements in HD 195987 and in other metal-poor stars 
cannot be ignored in the models, although their effect can be approximated in standard 
calculations by increasing the overall metallicity. Note, however, that this approach defeats 
the true purpose of the model comparisons, which is to test whether theory can reproduce 
all observed properties at the measured value of [Fe/H]. In the absence of a detailed chemical 
analysis for metal-poor stars, the constraint on models is considerably weakened because a 
satisfactory fit can usually be found by leaving the metallicity as a free parameter. The 
models by Bergbusch & VandenBerg (2001) incorporate the effect of a-enhancements, but 
regretably not the more refined boundary conditions used by Baraffe et al. (1998), which 
probably explains their failure to reproduce the My of the secondary and the effective tem- 
peratures of both stars. 

The example of HD 195987 emphasizes the importance of detailed spectroscopic analyses 
as a crucial ingredient for proper comparisons of metal-poor stars with stellar evolution 
theory. Not only is an accurate determination of [Fe/H] needed, but also abundances for 
other elements, particularly the a-elements that contribute the most to the opacities of these 
objects. To our knowledge the binary studied in this paper represents the most metal-poor 
object with such an analysis that has absolute masses determined to 2% or better. HD 195987 
is perhaps a favorable case for abundance determinations, however, because the secondary 
is faint enough that it does not interfere significantly with that analysis. This is seldom the 
case in double-lined eclipsing binaries that have accurately determined physical properties. 

While the tests described in previous sections appear to validate stellar evolution the- 
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ory, the constraint is actually not as strong as it could be because the absolute radii for the 
components are not directly measured (but see below). The stellar radius is a very sensi- 
tive measure of evolution (age), and the fact that it can be derived by purely geometrical 
means in binaries that are eclipsing is what makes such systems so valuable. From accurate 
observations of a number of systems that do eclipse, there is now clear evidence that while 
current models are quite successful in predicting the radiative properties of stars, they fail 
to reproduce the radii of stars in the lower main-sequence (see Popper 1997; Lastennet et al. 
1999; Clausen et al. 1999), predicting sizes that are too small. The discrepancy can be as 
large as 10-20%, as demonstrated recently for YY Gem AB and V818 Tau B by Torres & 
Ribas (2002). This may have a significant impact on ages inferred from these models. 

Incidentally, we note that the absolute radii that we derive in §5 for the components of 
HD 195987 based on the radiative properties of the system are also somewhat larger than 
predicted for stars of this mass by all of the models considered here, by roughly 10%. Though 
more uncertain and possibly not as reliable as the geometric determinations that might be 
obtained in eclipsing systems, our estimates are consistent with the trend described above. 

7. Final remarks 

New spectroscopic observations of the nearby moderately metal-poor double-lined spec- 
troscopic binary system HD 195987 ([Fe/H]~ —0.5, corrected for binarity), along with inter- 
feromctric observations that clearly resolve the components for the first time, have allowed us 
to obtain accurate masses {a < 2%), the orbital parallax, and component luminosities in V, 
H, and K. Both the detailed chemical composition (including the pattern of enhancement 
of the a-elements and heavier elements such as Ba and Eu) and the kinematics suggest that 
the system is a member of the thick disk population of our Galaxy, and an age of 10-12 Gyr 
or more (consistent with that conclusion) is inferred based on isochrone fits. 

The determination of these stellar properties places useful constraints on stellar evolu- 
tion theory for the lower main sequence. We have shown that none of the models considered 
here fits all the properties for both components simultaneously, although the discrepancies 
in each case can be understood in terms of the physical assumptions (convection prescrip- 
tions, boundary conditions, detailed chemical composition). Indications are that a model 
incorporating all the proper assumptions together (which is not done in any available series 
of calculations) would allow for a good match to the observations of HD 195987. While 
this may not be straightforward to do in practice due to the complexity of the problem, it 
would seem to be an obvious action item for theorists, which would enable observational as- 
tronomers to perform more useful tests for stars under 1 Mq. The ingredients identified here 
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as making the most significant difference in tfie fit for tfie mass regime of tfiis binary are: (a) 
tfie use of non-gray boundary conditions between the photosphere and the interior, based on 
modern model atmospheres that incorporate molecular opacity sources; (b) a mixing-length 
parameter close to that required for the Sun; and (c) the inclusion of the effect of enhanced 
abundances of the a-elements. 

Though it seems that all the radiative properties of HD 195987 can in principle be 
predicted accurately from current theory at the observed metallicity and ^-enhancement, 
there is no constraint in this particular case on the radius, a key diagnostic of evolution. 
Evidence cited earlier from other studies indicates that some adjustments are still needed in 
the models, which tend to underestimate the sizes of low-mass stars. 

Further progress on the observational side in testing models for metal-poor stars re- 
quires additional candidates suitable for accurate determinations of their physical properties. 
The recent paper by Goldberg ct al. (2002) presents a list of nearly 3 dozen double-lined 
spectroscopic binaries from a proper motion sample, many of which are metal-deficient. Im- 
provements in the sensitivity of ground-based interferometers in the coming years should 
allow several of them to be spatially resolved. Eclipsing binary candidates found in globular 
clusters in a number of recent studies are also valuable, but are fainter and the spectroscopy 
will be more challenging. In both cases detailed chemical analyses will be very important 
for a meaningful test of theory, as shown by the system reported here. 
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Fig. 1. — Relative orbit of HD 195987 on the plane of the sky, with the primary at the origin 
and the secondary shown at periastron. Motion is clockwise, and the dotted line indicates 
the line of nodes. The heavy line segments along the relative orbit indicate areas where 
we have phase coverage in our H- and i^-band PTI data (they are not separation vector 
estimates). Our data sample essentially all phases of the orbit well, leading to a reliable 
orbit determination. Component diameters are rendered to a scale three times larger than 
their actual size in relation to the orbit. 
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Fig. 2. — V"^ model for HD 195987. The top panel shows a direct comparison between four 
consecutive nights (24-27 June 1999) of calibrated K-hand V"^ data and the corresponding 
predictions computed from our "Full-fit" orbit model for HD 195987 (Table 6). The lower 
frame gives the residuals between the data and the model. The modified Julian Date is MJD 
= HJD - 0.5. 
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Fig. 3. — Radial velocity observations and spectroscopic orbit of HD 195987 from our 'Full- 
fit" solution (Table 6). The primary is represented with filled circles, and the center-of-mass 
velocity is indicated by the dotted line. The error bars for the primary are smaller than the 
point size. 
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Fig. 4. — Residuals from the orbital fit for HD 195987 (Table 6) as a function of orbital 
phase, for the H- and i^'-band visibilities {V^) and for the radial velocities, along with the 
corresponding histograms. 
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Fig. 5. — Comparison between the observations for HD 195987 and the Yale-Yonsei models 
by Yi et ah (2001). The left panels show isochrones for different ages, as labeled, and the 
right panels show the effect of a change in the metallicity. The infrared magnitudes from 
the original isochrones have been transformed to the CIT system (Bessell & Brett 1988). 
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Fig. 6. — The mass-luminosity diagram in the visual band, (a) Isochrones from the Yale- 
Yonsei models by Yi et al. (2001), as labeled, compared with the observations for HD 195987 
(circles) and other binary components in this mass regime (triangles). Infrared magnitudes 
have been converted to the CIT system, (b) Same as (a), but for the models by Baraffe et 
al. (1998), already on the CIT system. The agreement is much improved. 
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Fig. 7. — Comparison with the models by Baraffe et al. (1998) (left) and Siess, Forestini, & 
Dougados (1997) (right), showing the effect of age for solar composition, as in Figure 5. Also 
illustrated is the effect of a difference in the mixing-length parameter auL for the Baraffe et 
al. (1998) 12-Gyr isochrones (solid lines for uml = 1-0 and dot-dashed lines for uml = 1-9), 
and of a change in the overshooting parameter for the Siess, Forestini, & Dougados 
(1997) 12-Gyr isochrones (solid lines for ctov = 0.0 and dot-dashed lines for = 0.2ifp). 
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Fig. 8. — Effect of the light from the secondary in a binary, as a function of the magnitude 
difference in the visual band, (a) Changes in the color indices in the Johnson system, 
including the ultraviolet excess with and without correction for the guillotine (Sandage 1969); 
(b) Impact on the metallicity determination when using the ultraviolet excess 5{U—B) and the 
calibration by Marsakov & Shevelev (1988) (solid line), or the normalized excess 5{U —B)q,q 
(Sandage 1969) and calibrations by Carney (1979) (dashed lines). 
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Fig. 9. — Effect of the light from the secondary in a binary, as a function of the magnitude 
difference in the visual band, (a) Changes in the color indices in the Stromgren system. 
The indices represented with solid lines are used in the metallicity calibration by Schuster & 
Nissen (1989) adopted by Wyse & Gilmore (1995), and the ones shown as dashed lines are 
used in the calibration by Olsen (1984); (b) Impact on the metallicity determination when 
using the calibration by Schuster & Nissen (1989) (solid line) and by Olsen (1984) (dashed). 
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Fig. 10. — Effect of the light from the secondary as a function of the magnitude difference 
in the visual band, (a) Changes in the color indices, as labeled; (b) Bias in the temperature 
determinations using B — V, V — K, and b—y, and the calibrations by Carney et al. (1994). 
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Fig. 11. — Isochrone fits based on the models by Bergbusch & VandenBerg (2001) that 
include a-element enhancement, a) The fit shown is to the primary value of My only (see 
text), for an age of 12 Gyr and a fixed value of [a/Fe]= +0.36. The iron abundance from the 
fit is [Fe/H]= —0.45. A similarly good fit is obtained for an age of 10 Gyr and [Fe/H]= —0.52. 
Also represented for comparison is a model for the same [Fe/H] but with [a/Fe]= 0.0 (dashed 
line), b) Effective temperature predictions from the same models. 
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Table 1. Radial velocity measurements of HD 195987 and 

residuals from the final fit. 
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Table 1 — Continued 



HJD 


RVa 


RVb 


(O 


-C)a 


(O 


-C)b 


Orbital 


(2,400,000+) 


(kms^^) 


(kms^^) 


(km s ) 


(km s ^) 


Phase 


51380.600 


-25.37 


+20.08 




+0.18 




+0.91 


0.4673 


51413.642 


+28.58 


-49.05 




+0.00 




+0.49 


0.0437 


51424.724 


-10.20 


-0.90 




-0.19 




-0.34 


0.2371 


51445.726 


-24.67 


+18.91 




-0.26 




+1.18 


0.6035 


51464.487 


+21.74 


-41.97 




-0.10 




-0.99 


0.9308 


51466.675 


+29.24 


-51.29 




-0.10 




-0.79 


0.9689 


51468.667 


+32.45 


-52.45 




+0.51 




+1.35 


0.0037 


51476.569 


+7.09 


-21.63 




+0.06 




+0.54 


0.1415 


51521.513 


+20.77 


-42.93 




+0.15 




-3.50 


0.9256 


51525.510 


+32.26 


-53.03 




+0.48 




+0.57 


0.9953 


51530.440 


+20.79 


-42.20 




-0.28 




-2.20 


0.0813 


51618.884 


-23.54 


+18.30 




+0.06 




+1.60 


0.6243 


51640.847 


+31.85 


-53.88 




-0.06 




-0.12 


0.0074 


51667.839 


-25.55 


+18.42 




+0.16 




-0.95 


0.4783 


51682.802 


-14.67 


+10.21 




+0.44 




+4.29 


0.7393 


51699.727 


+29.26 


-51.44 




-0.63 




-0.24 


0.0346 


51712.696 


-12.89 


+4.55 




+0.13 




+1.28 


0.2609 


51788.531 


-24.41 


+18.05 




+0.59 




-0.43 


0.5838 


51808.604 


+22.50 


-42.91 




-0.10 




-0.97 


0.9340 


51815.609 


+26.25 


-47.50 




-0.16 




-0.73 


0.0562 


51842.553 


-26.18 


+17.75 




-0.29 




-1.85 


0.5263 


51866.503 


+24.55 


-46.12 




-0.28 




-1.34 


0.9441 


51874.448 


+20.28 


-39.25 




-0.47 




+0.35 


0.0827 


51921.441 


+14.54 


-33.52 




-0.39 




-1.31 


0.9025 


52004.892 


-21.96 


+13.64 




-0.47 




-0.38 


0.3583 



Table 2. Spectroscopic orbital solutions for HD 195987. 



Element This paper Imbert (1980) Latham et al. (1988) Duquennoy & Mayor (1991) Goldberg et al. (2002) 



Adjusted quantities 



P (days) 


57.32161 


± 


0.00034 


57.3210 


± 


0.0014 


57.325 


± 


0.009 


57.3240 


± 


0.0013 


57.3199 


± 


0.0065 


7 (km s^-'-) 


-5.867 


± 


0.038 


-5.630 


± 


0.075 


-5.56 


± 


0.12 


-6.13 


± 


0.07 


-5.59 


± 


0.12 


i^A (kms-i) 


28.944 


± 


0.046 


28.80 


± 


0.12 


28.89 


± 


0.20 


28.73 


± 


0.10 


29.16 


± 


0.19 


JCb (kms-i) 


36.73 


± 


0.21 




















34.6 


± 


1.2 


e 


0.3103 


± 


0.0018 


0.305 


± 


0.003 


0.316 


± 


0.006 


0.306 




0.003 


0.3083 


± 


0.0060 


t^A (dcg) 


357.03 


± 


0.35 


356.8 


± 


0.7 


355.9 


± 


0.8 


356.8 


± 


0.7 


357.1 


± 


1.1 


T (HJD-2,400,000)^. . 


49404.825 


it 


0.045 


49404.6 


± 


0.1 


49404.8 


± 


0.3 


49404.933 


± 


0.097 


49404.78 


± 


0.16 


lerived quantities 
































/(M) (M©) 








0.122 


± 


0.002 


0.1225 


± 


0.0026 


0.1219 


± 


0.0017 








MAsin^i {Mq) 


0.808 


± 


0.010 




















0.721 


± 


0.052 


Mesin^i (M©) 


0.6,369 


± 


0.0046 




















0.607 


± 


0.024 


q = Mb/Ma 


0.7881 


± 


0.0047 




















0.84 


± 


0.03 


OA sini (10^ km) 


21.689 


± 


0.036 


21.6 


± 


0.1 


21.65 


± 


0.16 


21.57 


± 


0.10 








Ob sin j (10^ km) 


27.52 


± 


0.16 


























asini (R©) 


70.70 




0.24 




















68.7 


± 


1.3 



Other quantities pertaining to the fit 

Afobs 73+73 60 31 64 28+28 

Time span (days) 6458 758 1003 4365 2653 

(TA (kms-i)!? 0.30 0.54 0.68 0.54 0.63 

(TB (kms-i)'? 1.38 ••• ••• ••• 4.51 



'^The time of periastron passage in the four solutions from the literature has been shifted by an integer number of cycles to the epoch derived in this paper, 
for comparison purposes. 

''Root mean square residual from the fit. 
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Table 3. Summary of the relevant parameters of the cahbration objects adopted for our 

analysis of the PTI visibihties. 







V 


K 




Adopted Model Diameter^ 


01)jeet Name 


S]")ectral Ty])e 


(mag) 


(mag) 


(dog) 


(mas) 


HD 195194 


G8 III 


7.0 


4.8 


2.7 


0.67 ± 0.10 


HD 200031 


G5 111 


6.8 


4.7 


6.0 


0.56 ± 0.05 


HD 177196 


A7 V 


5.0 


4.5 


17 


0.62 ± 0.10 


HD 185395 


F4 V 


4.5 


3.5 


17 


0.84 ± 0.08 



'^Angular separation from HD 195987. 

*^The angular diameters were determined from effective temperature and bolometric flux 
estimates based on archival broad-band photometry, and visibility measurements with PTI. 
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Table 4. Visibility measurements of HD 195987 in the H 
band, and residuals from the final fit. 



HJD 








u 


V 


Orbital 


(2,400,000+) 


y2 


(Jy2 


(0-C)v.2 


(10«A) 


(10«A) 


Phase 


51745.813 


0.322 


0.135 


-0.123 


-30.683 


-57.824 


0.8386 


51745.815 


0.407 


0.127 


-0.050 


-30.198 


-58.251 


0.8386 


51745.827 


0.639 


0.259 


+0.123 


-27.951 


-59.560 


0.8388 


51745.828 


0.608 


0.174 


+0.083 


-27.624 


-59.703 


0.8389 


51745.894 


1.051 


0.339 


+0.166 


-12.316 


-65.245 


0.8400 


51745.901 


0.920 


0.272 


+0.006 


-10.291 


-65.561 


0.8401 


51745.941 


1.025 


0.414 


+0.067 


+0.135 


-66.600 


0.8408 


52103.920 


0.743 


0.231 


+0.015 


-44.026 


-26.884 


0.0859 


52103.922 


0.740 


0.031 


-0.019 


-43.768 


-27.217 


0.0859 


52103.946 


1.054 


0.341 


+0.080 


-39.373 


-31.461 


0.0864 


52175.716 


0.584 


0.037 


+0.001 


-45.253 


-25.448 


0.3384 


52175.719 


0.638 


0.025 


-0.048 


-44.871 


-25.984 


0.3385 


52175.729 


0.911 


0.084 


-0.039 


-43.101 


-27.937 


0.3386 


52175.740 


0.956 


0.050 


+0.123 


-41.183 


-29.784 


0.3388 


52175.742 


0.846 


0.047 


+0.073 


-40.779 


-30.131 


0.3389 


52175.775 


0.481 


0.038 


-0.014 


-33.720 


-35.322 


0.3394 


52175.777 


0.499 


0.033 


—0.077 


—33.055 


—35.692 


0.3395 


52175.788 


0.854 


0.101 


+0.004 


-30.364 


-37.135 


0.3397 


52175.799 


0.905 


0.094 


+0.029 


-27.430 


-38.497 


0.3399 


52177.609 


0.980 


0.063 


+0.049 


-50.012 


-4.245 


0.3714 


52177.616 


0.977 


0.088 


+0.032 


-50.413 


-5.798 


0.3716 


52177.662 


0.789 


0.066 


+0.034 


-49.975 


-15.690 


0.3724 


52177.671 


0.977 


0.132 


+0.008 


-49.397 


-17.432 


0.3725 


52177.710 


0.534 


0.106 


+0.072 


-45.059 


-25.216 


0.3732 


52177.712 


0.662 


0.086 


+0.109 


-44.757 


-25.646 


0.3732 


52177.730 


0.944 


0.194 


+0.012 


-41.940 


-28.921 


0.3735 


52177.738 


0.651 


0.139 


-0.025 


-40.397 


-30.397 


0.3737 


52177.740 


0.643 


0.089 


+0.045 


-39.995 


-30.719 


0.3737 


52179.617 


0.972 


0.059 


+0.028 


-50.720 


-7.199 


0.4065 


52179.625 


0.978 


0.061 


+0.075 


-50.815 


-8.877 


0.4066 


52179.662 


0.692 


0.072 


+0.045 


-49.836 


-16.701 


0.4072 


52179.669 


0.942 


0.082 


+0.019 


-49.322 


-18.258 


0.4074 


52179.704 


0.373 


0.059 


+0.027 


-45.322 


-25.308 


0.4080 


52179.706 


0.446 


0.027 


+0.036 


-45.064 


-25.641 


0.4080 


52179.713 


0.802 


0.077 


+0.061 


-43.940 


-27.019 


0.4081 


52179.731 


0.736 


0.068 


-0.029 


-40.884 


-30.179 


0.4085 


52179.739 


0.415 


0.058 


-0.026 


-39.393 


-31.485 


0.4086 
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Table 5. Visibility measurements of HD 195987 in the K 
band, and residuals from the final fit. 



HJD 














u 




V 




Orbital 


(2,400,000+) 




<Ty2 


(0-C)^2 


(106 A) 


(lO^A) 
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.3342 
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Table 5 — Continued 



HJD 








It 


V 


Orbital 


(2,400,000+) 


y2 


(J 2 




(10^ A) 


(10^ A) 


Phase 


51716.953 


1.015 


0.075 


+0.135 


-12.895 


-47.520 


0.3351 


51716.963 


0.935 


0.101 


+0.078 


-11.014 


-48.043 


0.3353 


51716.990 


0.565 


0.059 


-0.052 


-5.850 


-48.819 


0.3358 


51717.001 


0.570 


0.065 


+0.065 


-3.743 


-49.051 


0.3359 


51723.876 


0.874 


0.067 


+0.025 


-22.044 


-42.796 


0.4559 


51723.890 


0.873 


0.075 


+0.044 


-20.086 


-43.970 


0.4561 


51723.919 


0.284 


0.046 


-0.044 


-15.387 


-46.184 


0.4566 


51723.933 


0.426 


0.061 


+0.024 


-12.929 


-46.946 


0.4569 


51723.949 


0.715 


0.066 


+0.045 


-9.883 


-47.728 


0.4572 


51723.963 


0.866 


0.080 


+0.055 


-7.236 


-48.294 


0.4574 


51723.974 


0.822 


0.224 


+0.014 


-5.106 


-48.406 


0.4576 


51723.988 


0.686 


0.089 


+0.001 


-2.333 


-48.729 


0.4579 


51723.996 


0.710 


0.227 


+0.105 


-0.918 


-48.809 


0.4580 


51741.797 


0.976 


0.055 


+0.048 


-36.915 


-2.364 


0.7685 


51741.804 


0.735 


0.090 


-0.153 


-37.236 


-3.540 


0.7687 


51741.902 


0.397 


0.019 


-0.024 


-33.607 


-18.647 


0.7704 


51741.909 


0.559 


0.036 


+0.019 


-32.877 


-19.660 


0.7705 


51741.944 


1.082 


0.038 


+0.094 


-28.023 


-24.126 


0.7711 


51741.951 


1.042 


0.041 


+0.066 


-26.960 


-24.913 


0.7712 


51744.796 


0.254 


0.024 


-0.011 


-37.043 


-3.596 


0.8209 


51744.803 


0.235 


0.032 


-0.021 


-37.251 


-4.670 


0.8210 


51744.851 


0.200 


0.035 


-0.003 


-36.705 


-12.146 


0.8218 


51744.858 


0.191 


0.020 


-0.013 


-36.328 


-13.234 


0.8219 


51744.928 


0.484 


0.047 


-0.024 


-29.065 


-23.065 


0.8232 


51744.935 


0.542 


0.077 


-0.023 


-28.034 


-23.957 


0.8233 


51746.822 


0.206 


0.045 


+0.000 


-21.122 


-44.293 


0.8562 


51746.824 


0.226 


0.049 


+0.019 


-20.782 


-44.390 


0.8562 


51746.836 


0.253 


0.044 


+0.042 


-19.017 


-45.391 


0.8564 


51746.839 


0.287 


0.086 


+0.076 


-18.522 


-45.707 


0.8565 


51746.893 


0.301 


0.023 


+0.047 


-8.869 


-48.761 


0.8574 


51746.895 


0.300 


0.099 


+0.043 


-8.427 


-48.856 


0.8575 


51752.793 


0.460 


0.032 


-0.008 


-37.534 


-6.427 


0.9604 


51752.799 


0.596 


0.040 


+0.017 


-37.575 


-7.511 


0.9605 


51752.806 


0.701 


0.058 


+0.007 


-37.545 


-8.587 


0.9606 


51752.829 


0.955 


0.075 


-0.015 


-36.897 


-12.119 


0.9610 


51752.835 


0.999 


0.082 


+0.011 


-36.571 


-13.156 


0.9611 


51752.842 


0.944 


0.066 


-0.024 


-36.237 


-14.146 


0.9612 


51752.872 


0.468 


0.048 


-0.040 


-33.585 


-18.607 


0.9617 


51791.675 


0.927 


0.084 


-0.007 


-24.303 


-41.631 


0.6387 


51791.683 


0.925 


0.068 


+0.082 


-23.256 


-42.516 


0.6388 


51794.671 


0.282 


0.050 


+0.019 


-23.622 


-41.847 


0.6909 


51794.714 


0.995 


0.096 


+0.064 


-17.173 


-45.574 


0.6917 


51794.780 


0.420 


0.051 


+0.047 


-5.164 


-48.662 


0.6928 


51831.634 


0.981 


0.282 


+0.117 


-13.669 


-47.076 


0.3358 


51831.644 


0.957 


0.364 


+0.077 


-11.832 


-47.527 


0.3359 


51832.645 


0.618 


0.038 


-0.016 


-11.295 


-48.013 


0.3534 


51832.655 


0.513 


0.057 


+0.026 


-9.402 


-48.524 


0.3536 


51832.687 


0.327 


0.032 


-0.006 


-3.219 


-49.407 


0.3541 
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Tablc 5 — Continued 



HJD 








u 


V 


Orbital 


(2,400,000+) 


y2 


(7y'2 


{0-G)y1 


(106 A) 


(lO^A) 


Phase 


51832.703 


0.388 


0.051 


-0.017 


+0.028 


-49.456 


0.3544 


51832.741 


0.613 


0.079 


-0.006 


+7.522 


-48.791 


0.3551 


51832.749 


0.747 


0.185 


+0.080 


+9.092 


-48.683 


0.3552 


51833.587 


0.779 


0.073 


+0.054 


-20.747 


-44.399 


0.3698 


51833.594 


0.797 


0.076 


-0.030 


-19.686 


-44.961 


0.3700 


51833.628 


0.542 


0.064 


-0.062 


-13.791 


-47.277 


0.3706 


51833.636 


0.370 


0.090 


-0.100 




-47.827 


0.3707 


51833.671 


0.404 


0.065 


+0.020 


-5.797 


-49.071 


0.3713 


51833.678 


0.456 


0.055 


+0.005 


-4.458 


-49.294 


0.3714 


51833.713 


0.875 


0.229 


+0.136 


+2.534 


-49.312 


0.3720 


51833.720 


0.959 


0.244 


+0.184 


+3.949 


-49.290 


0.3722 


51850.576 


0.299 


0.206 


+0.092 


-35.890 


-14.261 


0.6662 


51850.587 


0.333 


0.143 


-0.066 


-35.123 


-16.014 


0.6664 


51850.589 


0.774 


0.457 


+0.331 


-34.953 


-16.216 


0.6664 


51850.627 


0.842 


0.286 


+0.130 


-30.666 


-21.528 


0.6671 


51850.636 


0.425 


0.144 


+0.020 


-29.413 


-22.756 


0.6673 


51851.630 


0.586 


0.128 


-0.095 


-29.867 


-22.233 


0.6846 


51851.632 


0.588 


0.039 


-0.025 


-29.662 


-22.568 


0.6846 


52120.799 


0.244 


0.015 


-0.020 


-37.790 


-8.720 


0.3804 


52120.806 


0.435 


0.028 


+0.017 


—37.598 


—9.746 


0.3805 


52120.821 


0.933 


0.039 


+0.063 


-37.127 


-12.097 


0.3807 


52120.827 


1.037 


0.034 


+0.058 


-36.800 


-13.153 


0.3809 


52120.850 


0.517 


0.030 


-0.009 


-35.175 


-16.534 


0.3812 


52154.730 


0.700 


0.112 


+0.016 


-17.204 


-46.179 


0.9723 


52154.753 


1.056 


0.075 


+0.153 


-13.188 


-47.569 


0.9727 


52154.776 


1.025 


0.099 


+0.114 


-8.897 


-48.621 


0.9731 


52156.711 


0.447 


0.017 


+0.023 


-19.236 


-44.927 


0.0069 


52156.741 


0.924 


0.032 


+0.086 


-14.316 


-46.944 


0.0074 


52156.767 


0.949 


0.025 


+0.030 


-9.421 


-48.279 


0.0078 


52156.788 


0.715 


0.025 


-0.046 


-5.397 


-48.945 


0.0082 


52180.651 


0.373 


0.026 


-0.017 


-37.374 


-11.148 


0.4245 


52180.658 


0.657 


0.027 


+0.021 


-37.059 


-12.252 


0.4246 


52180.671 


1.034 


0.036 


+0.052 


-36.298 


-14.354 


0.4249 


52180.678 


0.955 


0.031 


+0.013 


-35.826 


-15.419 


0.4250 


52180.692 


0.458 


0.034 


-0.028 


-34.620 


-17.493 


0.4252 


52180.699 


0.280 


0.018 


-0.007 


-33.960 


-18.397 


0.4253 


52180.713 


0.310 


0.017 


+0.004 


-32.398 


-20.408 


0.4256 


52180.720 


0.543 


0.021 


+0.018 


-31.499 


-21.269 


0.4257 



Table 6. Orbital parameters for HD 195987. 



Parameter HIPPARCOS^ 

Adjusted quantities 

P (days) 57.3240 

a (mas)l? 5.24 ± 0.66 

7 (km s"-*-) 

Ka (kms-i) 

Kb (kms-i) 

e 0.3060 

i (deg) 89.50 ± 8.43 

LOA (deg) 356.8 

f2 (deg)<; 327.66 ± 7.56 

T (HJD-2,400,000) .... 43327.589^1 

Aii'ciT (mag) 

A-ffciT (mag) 

Derived quantities 

TT (mas)"r 45.30 ± 0.46 

Other quantities pertaining to the fit 

iVobs (RV) 

iVobs (^2) 

ay2(H/Ky. 

(TRV (A/B, kms-i)f. ... 



y2 Only 



RV Only 



Full Fit 



57.3298 


± 


0. 


.0035 


57.32161 


± 


0.00034 


57.32178 


± 


0.00029 


15.368 


± 


0, 


.028 








15.378 


± 


0.027 










-5.867 


± 


0.038 


-5.841 




0.037 










28.944 


± 


0.046 


28.929 


± 


0.046 










36.73 


± 


0.21 


36.72 




0.21 


0.30740 


± 


0, 


.00067 


0.3103 


± 


0.0018 


0.30626 


± 


0.00057 


99.379 


± 


0, 


.088 








99.364 


± 


0.080 


358.89 


± 


0, 


.53 


357.03 


± 


0.35 


357.40 


± 


0.29 


335.061 


± 


0, 


.082 








334.960 


± 


0.070 


51354.000 


± 


0, 


.069 


49404.825 


± 


0.045 


51353.813 




0.038 


1.063 


± 


0, 


.031 








1.056 


± 


0.013 


1.18 


± 


0, 


.16 








1.154 


± 


0.065 
















46.08 


± 


0.27 



37+134 
0.0599 / 0.0613 



73+73 



0.30 / 1.38 



73+73 
37+134 
0.0592 / 0.0617 
0.31 / 1.41 



''The elements P, e, loa, and T in this solution were adopted from Duquennoy & Mayor (1991) and held fixed. 

''For the HIPPARCOS solution this is the semimajor axis of the photocenter (a) rather than the relative semimajor axis. 

"^Position angle of the ascending node, i.e., the node at which the secondary is receding. 

"^This value listed in the HIPPARCOS catalog is erroneous; it should be 43328.589 (cf Duquennoy & Mayor 1991). 

''The determination from the HIPPARCOS data is our revised value (see §4), and the full- fit determination is the orbital parallax. 

^Root mean square residual from the fit. 
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Table 7. Physical parameters of HD 195987. 



Parameter Primary Secondary 



Mass (M0) 0.844 ± 0.018 0.6650 ± 0.0079 

Teff (K) 5200 ± 100 4200 ± 200 

TTorb (mas) 46.08 ± 0.27 

Dist (pc) 21.70 ± 0.13 

My (mag) 5.511 ± 0.028 7.91 ± 0.19 

Mh"" (mag) 3.679 ± 0.037 4.835 ± 0.059 

Mk"" (mag) 3.646 ± 0.033 4.702 ± 0.034 

V-K^ (mag) 1.865 ± 0.039 3.21 ± 0.19 



^In the CIT system. 



